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c^ ■ Abstract 



A precise test of the theory of stehar evolution can be performed by mea- 
^ . suring the average difference in energy between the neutrino hue produced 

CN ■ by ^Be electron capture in the solar interior and the corresponding neutrino 



o 



1 



line produced in a terrestrial laboratory. This energy shift is calculated to be 



o 

^\ . 1.29 keV (to an accuracy of a few percent) for the dominant ground-state to 

Oj. ground-state transition. The energy shift is approximately equal to the aver- 

o 

d ■ over the solar interior with a weighting factor equal to the locally-produced 



age temperature of the solar core, computed by integrating the temperature 



^Be neutrino emission. Therefore, a measurement of the energy shift is a 
measurement of the central temperature distribution of the sun. 

The energy profile of the '^Be line is derived analytically and is evaluated 
numerically. The line shape is asymmetric: on the low-energy side, the line 
shape is Gaussian with a half-width at half-maximum of 0.6 keV and on 
the high-energy side, the line shape is exponential with a half-width at half- 
maximum of 1.1 keV. The effective temperature of the high-energy exponential 
tail is 15 X 10^ K. The energy profile of the ^Be neutrino line should be 
taken into account in calculations of vacuum neutrino oscillations and of the 
absorption cross section for ^Be solar neutrinos incident on '^Li nuclei. 



The characteristic modulation of the Be hne shape that would be caused 
by either vacuum neutrino oscillations or by matter-enhanced (MSW) neu- 
trino oscillations is shown to be small. Other frequently-discussed weak inter- 
action solutions to the solar neutrino problem are also not expected to change 
significantly the line profile. 
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I. INTRODUCTION 



A. Background 



The original motivation |[l|] for performing solar neutrino experiments was to learn more 
about how the sun shines. When the first observational results from the chlorine experi- 
ment became available [Q, the focus shifted [^ from learning about the interior of the sun 
to trying to determine if the discrepancy between calculation and observation was due to 
inadequacies in the astrophysics or to new weak interaction physics. Directed towards this 
goal, four experiments are being performed [^|-0] to measure the fluxes of neutrinos that are 
produced by nuclear fusion reactions in the solar interior. An additional four solar neutrino 
experiments P-|TT[] — designed to determine if new physics or new astronomy is required — are 
being developed. 

Progress toward the goal of testing solar theory by measuring solar neutrinos has been 
complicated by what may be the discovery of new weak interaction physics. A comparison 
of two of the existing experiments, the chlorine [HI and the Kamiokande H] experiments. 



suggests [|T2[ the existence of a physical process that changes — in a way that depends upon 
neutrino energy — the fraction of electron-type neutrinos that reach detectors on earth after 
being created in the core of the sun. 

At first glance, the terrestrially-observed flux of ^B neutrinos is encouragingly close to 
the flux calculated on the basis of the standard solar model and the standard electroweak 
theory, especially considering the sensitivity of the predicted flux to details of the stellar 
physics. However, the discrepancy between calculation and observation is significant since 
the theoretical uncertainties are smaller than the difference between what is measured and 



what is predicted ||T3Hl5|. This discrepancy has stimulated a number of imaginative and at- 



tractive possible interpretations in terms of new physical processes [0-22]. Until the effects 
of these proposed new physical processes are either established or rejected experimentally, 
quantitative astronomical inferences from the measured magnitudes of the solar neutrino 



fluxes will be limited. 

B. The ^Be Energy Profile 

The purpose of this article is to draw attention to precise new predictions of the theory of 
stellar evolution that can be tested with the aid of future solar neutrino experiments. These 
new predictions are based upon the calculated shape of the energy spectrum of neutrinos 
produced by ^Be electron capture in the solar interior. The reaction in question is 

e" +^ Be -^ ^Li + z/^ . (1) 

In the interior of the sun, most of the electrons are captured from continuum (unbound) states 



23| , |2^ . The electron-capture reaction shown in ([^) produces a neutrino line because the ^Li 



nucleus in the final state has a mass that is much greater than the energy of the neutrino. 
The recoiling nucleus takes up a significant amount of momentum, but only a negligible 
amount of energy. The focus of this paper is on the broadening by thermal effects of the 
neutrino line produced in the sun. 

This paper shows that the line shape for Eq. (|l|) reflects accurately the temperature 
distribution in the interior of the sun. As we shall see, the difference between the average 
energy of the neutrino line emitted in the sun and in the laboratory is approximately equal to 
the central temperature of the sun (see especially Sec. [VB| and Sec. |VIB| ). The predicted 1.3 
keV increase in the average energy of the solar ^Be neutrino line relative to laboratory decays 
is the simplest quantity to measure that directly reflects the solar temperature distribution. 

In addition to a shift in the average neutrino energy, the shape of the line profile for ^Be 
electron capture in the sun is different from the line profile that would be observed for a 
laboratory source of ^Be neutrinos. The shift in average energy and the change in the shape 
of the line profile are both caused by the high temperatures in the core of the sun where the 
neutrinos are produced. The high solar temperatures produce significant thermal energies 
for continuum electrons and their capturing nuclei, Doppler shifts of the emitting nuclei, a 



high degree of ionization of solar ^Be ions, and a difference in atomic binding energies for 
solar and laboratory ^Be atonis.Q 

Figure 1 shows the two nuclear transitions that occur when ^Be captures an electron in 
the laboratory. The neutrino energy corresponding to the transition from the ground-state 
of ^Be to the ground-state of ^Li will be denoted by gLab(g-S-)) the energy corresponding to 
the transition to the first excited state of ^Li will be denoted by gLab(ex.s.). In what follows, 
I will discuss both transitions on an equal basis since the physical processes determining 
the neutrino line shape are the same in both cases. However, the transition to the ground 
state of '^Li is more easily studied experimentally because it has a higher energy and a larger 
branching ratio p^. Both the higher energy and the larger branching ratio of the ground- 
state to ground-state transition contribute to making the expected rate for this transition 
faster by an order of magnitude than for the ground-state to excited-state transition. 

For convenience, I will refer to the ground-state to ground-state transition as "ground- 
state" capture and will refer to the ground-state to excited-state transition as "excited- 
state" capture. In an additional effort to avoid cumbersome phrasing, I will often refer in 
the singular to "the" line profile or to "the" energy shift, when I mean the line profile or 
the energy shift for both ground-state capture and excited-state capture. 

Before proceeding to the calculations, it is helpful to think about the following question. 
Why is the effect of the solar environment on the shape of the energy spectra for continuum 
beta-decays (one particle in the initial state, three particles in the final state, e.g., ^B decay) 
different from the effect of solar conditions on the profile shape for the two-body electron 
capture reactions? This question has a simple physical answer. 

If the simplest version of the standard electroweak model is correct, then the shape of 



^ After the initial version of the present work was accepted for publication |2^, my attention was 
directed by J. Pantaleone to two important papers [^. In these two papers, the effect of the line 
width of the ^Be and the pep lines on neutrino oscillations was calculated. 



the energy spectrum for electron-type neutrinos from continuum beta-decays in the sun, 
such as the '^B beta-decay, is independent of the conditions in the sun to an accuracy of 
1 part in 10^ |^. The invariance is a result of the fact that, in the center-of-momentum 
frame, a decaying ®B nucleus has no kinetic energy. In the laboratory frame, terms of order 
the velocity of the nucleus, f (^B)/c, cancel out because there are as many ^B nuclei moving 
toward the observer as there are moving away from the observer. The first-order Doppler 
effects vanish because at each point within the broad ^B continuum the neutrino energy is 
spread out symmetrically by a small amount. The largest potentially-observable effects of 
the solar temperature on the observed ^B neutrino energy spectrum are only second order 
terms, ~ f (^B)^/c^ ||2^. This implication of standard theory will be tested by experiments 



|8|j9|jTl|| that will measure the shape of the *B neutrino energy spectrum and will compare 
the observed shape with the spectrum determined from terrestrial measurements. 

By contrast, in the two-body electron-capture reactions, the profile of a narrow neutrino 
emission line is broadened asymmetrically by the solar temperature. On the low-energy 
side of the neutrino line profile, the dominant effect is first-order Doppler (or Gaussian) 
broadening caused by the motion of the decaying nuclei. On the high-energy side of the 
observed profile, the dominant effect is an exponential broadening resulting from the center- 
of-momentum kinetic energy of the electron and the decaying nucleus (for a physical discus- 



sion of the profile shape see especially Sec. VI). This asymmetric broadening causes a shift 



of the average energy of the neutrino line as well as change in the shape of the line profile. 

C. Experimental Possibilities 

A number of experiments 0,^ have been proposed that would measure predominantly 
the z/g flux from ^Be electron capture in the sun, using detectors that are based upon 
neutrino-electron scattering. At the present time, the BOREXINO experiment ||T0| is the 
most advanced of these proposals and can, if recent estimates of the expected backgrounds 
are correct, measure the flux of ^Be neutrinos. Radiochemical detectors of the flux of ^Be 



neutrinos have also been discussed |^, but these detectors do not give specific information 



on neutrino energies. The radiochemical detectors are efficient for measuring the total ffux 
of electron-type neutrinos above a specified energy threshold, but are not useful for studying 
the thermal effects investigated in this paper. Most recently, the use of a high resolution 
LiF bolometric detector of ^Be neutrinos has been discussed 



Detectors have been developed |3^ for a variety of other applications, including dark 
matter searches, the observation of double beta decay, and x-ray astronomy, that have the 
energy resolution and the sensitivity that are required to study the energy spectrum of the 
^Be neutrino line. The best-available detectors have energy resolutions, AE/E, of much 
better than 1%, but they are smaller than would be required for a full-scale solar neutrino 
detector. 

The most direct way to study the ^Be energy profile may be to detect neutrino absorption 
by nuclei, a process which leaves an electron and a recoiling nucleus in the final state. In 
these charged-current transitions, nearly all of the initial neutrino energy is transferred to 
the final-state electron (the nuclear recoil energy being small). In the neutrino-electron 
scattering experiments that are currently under development, the 1 keV width of the ^Be 
line is spread out over several hundred keV of electron recoil energy, since the neutrino and 
the electron share the final state state energy (see Figure 8.5 of [ll3|). 

In order to measure the predicted 1.3 keV (0.15%) energy shift via neutrino absorption, an 
energy resolution of order 1% to 0.1% is desirable, depending somewhat upon the absorption 
threshold. Consider, for specificity, a conceivable cryogenic experiment [Q that might 
be performed on ^^Br with an energy resolution of 1% and with a total of 10^ measured 
neutrino events. The energy released to the recoil electron would be about 400 keV (the 
reaction threshold is about 450 keV), so the average neutrino energy would be measured to 
an accuracy of about 0.1 keV. The experimental parameters assumed in the above discussion 
would permit one to measure the central temperature of the sun to an accuracy of about 
10%. In addition, a proposed high-pressure hehum gas detector [Q might well have sufficient 
energy resolution to measure the predicted energy shift. 



The requirements for a practical experiment may be achievable since the solar neutrino 
detectors currently under development are designed to detect several thousand events per 
year (albeit with much poorer energy resolution, typically ~ 10%). It would be valuable to 
calibrate the solar results by studying an intense laboratory source of ^Be neutrinos with 
the same detector as used in the solar observations. The work described in this paper was 
undertaken in the hope that it would stimulate an experiment that would measure the energy 
shift and perhaps other characteristics of the ^Be line profile, in somewhat the same way 
that the initial theory and the experimental results on solar neutrinos developed together 

i- 

D. Organization and Previous Work 

This paper is organized as described below. Section |I| presents calculations of the average 
neutrino energy release in the rest frame of the decaying particle when a ^Be nucleus captures 
an electron under laboratory conditions (see Sec. [11 A| ) and in a solar environment (in which 
most of the electrons are captured from continuum orbits, see Sec. [11 Bp . Section ^ describes 
the central calculation of this paper, an evaluation of the energy profile for the neutrino line 
emitted when ^Be nuclei capture continuum electrons that have a specified temperature. 
The following section. Sec. |I^, outlines the calculation of the energy profile for the small 
but significant fraction of the solar electron captures that occur from bound orbits. The 
numerical characteristics of the line shape, including the shift in average neutrino energy 
between laboratory and solar decays, the full-width at half-maximum of the line profile, and 
the lower-order moments of the energy distribution, are the subject of Sec. 0. The line profile 
is averaged, in this section, over the physical characteristics, including the temperature 
distribution, of detailed solar models. The numerical results are presented in Tables I and II 
and displayed in Figures 2 and 3. Section |VT| provides approximate analytic derivations of the 
low-energy half-width (0.56 keV) and the high-energy half-width (1.07 keV) of the neutrino 
energy profile and isolates the separate physical origins of these two features. This section 
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also presents a derivation of the shift in average neutrino energy between solar and laboratory 
decays that explains why the shift is the same for ground-state and for excited-state captures. 
Section |V11| evaluates the effects on the neutrino line profile of electrostatic screening, of 
gravitational redshifts, and of collisional broadening, and shows that these effects are small. 
The influence of possible new weak interaction physics, exemplified by vacuum oscillations 
and by the MSW effect, is the subject of Sec. |VIII| . Neither type of oscillation would affect 
significantly the measurable characteristics of the ^Be line profile. The energy shift and the 
asymmetric profile of the ^Be line do change the computed absorption cross section for ^Be 
solar neutrinos incident on a ^Li detector, as is shown in Section |I^. Section ^ summarizes 
the principal results of this paper. 

There have been a number of previous calculations of the total rate for ^Be electron 
capture in the sun over the past three decades p3|,p4|,p5|-|39|| . The present work on the ^Be 



line shape does not change the results of these prior calculations of the total capture rate, the 



standard value in general use still being the one given in reference [0. The only earlier work 
on the broadening of the ^Be neutrino line with which I am familiar concentrated either on 
the fraction of electron-capture neutrinos that were above threshold for the i/g + ^Li -^ ^Be+e 



reaction, taking account of center-of-momentum energies but not Doppler shifts B^Jc 



or 



on a crude estimate of the total width, based only on the Doppler shifts of the ^Be ions [|T3 



II. LABORATORY AND SOLAR ENERGIES 

The average neutrino energy that is released when ^Be captures an electron in a terrestrial 
laboratory is different from the average energy that is released when ^Be captures an electron 
in the core of the sun. In a laboratory experiment, all of the captured electrons are initially 
bound in a ^Be atom, whereas in the solar interior most of the electron captures occur from 
continuum orbits |2^. Part of the difference in the average neutrino energy release, the part 
on which this section concentrates, is due, therefore, to the different atomic binding energies 
in the laboratory and in the sun. 
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Energy conservation implies, for electron captures either in the sun or in the laboratory, 
that 

^initial - ^finai = AM + K{e) + ^^(^Be) - KCU) - q + aCBe) - a(^Li) = 0, (2) 

where -^initial and -Efinai are the total energies in the initial and final states, K{e) and i^(^Be) 
are kinetic energies of the electron and of the ^Be nucleus in the initial state, K(^Li) is the 
kinetic energy of the recoiling ^Li nucleus, q is the neutrino energy, and a(^Be) and a(^Li) 
are the atomic binding energies. The quantity AM in Eq. (0) is the difference in rest mass 



energies [g^]; excluding atomic binding energies, 

AM = me + M(^Be) - M(^Li) = 862.10 keV. (3) 

For electron captures that occur in the laboratory, the initial kinetic energies, K{e) and 
i^(^Be), are both zero. However, these terms contribute numerically the largest amount to 
the difference between laboratory and solar neutrino release and will be calculated in Sees. 



Ill and IV 



In the following two subsections, I evaluate the neutrino energy using Eq. (|]) for labo- 
ratory decays (in Sec. |1I A| ) and for solar decays (in Sec. |1IB|) . I discuss in Sec. [VII A| the 



electrostatic energy of the screening cloud that surrounds the decaying nucleus. 

A. Laboratory Decays 

The laboratory neutrino energy, q^nh, satisfies the equation 

?Lab + qLj2MCU) = AM + aLab('Be) - aLab('Li), (4) 

where the difference in atomic binding energies is 

aLab('Be) - aLab('Li) = -0.195 keV. (5) 

Since the energy difference between an initial and a final state is independent of the choice 
of steps used to reach the final state, the difference in atomic binding energies is equal to 

10 



the difference between the sum of the successive ionization potentials of ^Be and the sum of 
the successive ionization potentials of ^Li. The difference in atomic binding energies given 
in Eq. (^ was calculated by subtracting the sum of the three ionization potentials of the 
lithium atom from the sum of the four ionization potentials of atomic beryllium, using the 
measured values for the ionization potentials |^ . 



Combining Eq. (H) and Eq. (H), one obtains 

me + M(^Be) - M(^Li) + aLab(^Be) - aLab(^Li) = 861.90 keV, (6) 



which is the tabulated mass difference of the neutral atoms [^. The small contribution of 
the ^Li recoil energy is 

g^(g.s.)/2M(^Li) = 0.057 keV ; q\ex.s.)/2MCU) = 0.011 keV, (7) 

for ground-state and excited state transitions, respectively. Inserting Eq. (H) and Eq. (0) in 
Eq. (^ , we find the neutrino energies for the laboratory transitions, 

gLab(g.s.) = 861.84 keV, (8) 

and 

gLab(ex.s.) = 384.28 keV, (9) 

that are shown in Fig. 1. 

B. Solar Decays 

This subsection begins with a calculation of the neutrino energy that is released in the 
limiting case in which a continuum electron with zero kinetic energy is captured by a ^Be 
nucleus at rest. I then calculate the neutrino energy that is emitted when a bound electron is 
captured from a stationary ^Be nucleus. I evaluate the influence of non-zero kinetic energies 
of the electron and of the ^Be ion in Sees. |T| and ^V\. 
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For the special case of continuum electron capture at rest, the neutrino energy release is 
the same as for laboratory capture except that the nucleus is assumed stripped of all bound 
electrons. Referring to Eq. (^, one can write: 

gcont.star = Qhah " KabC^Be) - ai^^iJU)]. (10) 

Therefore, 

gcont,star = AM - gLt,star/2M(^Li), (11) 

where the difference, AM, in rest mass energies is given by Eq. (^. Thus 

gcont,star(g.S.) = 862.04 kcV . (12) 

I next consider the average neutrino energy that is released, (jbound.star , when an electron 
is captured from a bound orbit. The only difference between gbound,star and the previously 
calculated gcont.star is the atomic binding energy of the electrons. Therefore, one can write 

Q'bound.star = Q'cont.star + (C^star( Be) — Ostari ■L'l)/) (13) 

where the angular brackets denote an average over the sun. The K-shell binding energy has 



been determined for solar conditions by Iben and his collaborators pj] using a variational 
calculation. For the case in which one electron is bound to the ^Be nucleus, the result can be 
expressed in terms of the relative K-shell binding, a^, in the sun compared to the laboratory 
value, where a^ is defined by the following equation: 

aK,star('Be) = -2lQ.QaR eV. (14) 

At the peak of the ^Be solar neutrino emission, a^ = 0.25 and aK,star(^Be) = —0.06 keV 
(cf. discussion in Sec. |V|). This binding energy is sufficiently small that it is not necessary, 
for our purposes, to calculate aK,star(^Be) to high precision. However, an accurate calculation 
was carried out without difficulty and the results of this calculation are described below. The 
average binding energies were obtained for two standard solar models [jl5| (see also Sec. [VA| ), 
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one that included, and one that did not include, helium diffusion. Using the tabulated results 
of the variational-principle calculation [^, a convenient interpolation formula was derived 



for an in terms of the local density, temperature, and chemical composition. Weighting 
the value of aR calculated at a given radial distance from the center of the sun by the ^Be 
neutrino flux at that radius, the average relative K-shell binding energy for the standard 



solar model ||15[ including helium diffusion is 

an ~ 0.255, (15) 

and is 0.256 for the solar model that does not include helium diffusion. Thus, the average 
binding energy of a K-shell electron in the sun is 

(aK,star('Be)) = -0.055 keV. (16) 

The difference in atomic binding energies between two K-shell electrons in a beryllium atom 
in the sun and a single K-shell electron in a lithium atom in the sun may be estimated to 
more-than-adequate accuracy by scaling the K-shell energy given in Eq. (0) by the ratio 
of the square of the lithium and the beryllium nuclear charges. One can also estimate 
the energy difference by scaling (see l2^), using the same value of gr calculated above. 



the laboratory energy that is required to remove one electron from previously twice-ionized 
beryllium, leaving behind one electron bound to a lithium atom. Both methods give the 
same answer to the accuracy of interest here. One finds 

(a2K,star('Be) - aK,star('Li)) ~ -0.079 kcV. (17) 

One must average the results given in Eq. (^) and Eq. (|1^) over the stellar model taking 
account of the variable fraction of the decays that occur from ^Be atoms with one or with 
two bound electrons. Let pi and p2 be the respective probabilities that "^Be has one or two 
bound electrons that can decay by K-capture. Then the appropriate average can be written 



as 



/ (It:, \ (It -An /Pl«K,star C^Bc) + 2^2 (a2K,star (^Bc) - a2K,star (^Li)) 

(a ( Bej - a ( Lij),tar = ( '^^^^^ ^ ^^^^ 
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Convenient expressions for pi and p2 have been given by Iben and his collaborators p4 . 
Carrying out the average using the results for the relative ionization states obtained from 
the variational-principle calculation, I find for the atomic binding energy in the sun: 

(a(^Be) - a(^Li))star = -0.06 keV. (19) 

Inserting Eq. (|19D in Eq. (|l|), 

gbound,star(g.S.) = 861.98 keV. (20) 

By an analogous procedure, one can obtain the following results for the excited-state 
transition: 

gcont,star(ex.s.) = 384.43 keV, (21) 

and 

gbound,star(ex.s.) = 384.38 keV. (22) 

III. CAPTURE FROM CONTINUUM ORBITS 

This section presents a calculation of the energy profile for the neutrino line emitted 
when electrons are captured from continuum orbits by ^Be ions in the core of the sun. The 
energy profile is the probability that a ^Be ion captures an electron from the continuum and 
emits a neutrino of any specified energy. As noted earlier, continuum electron capture is the 
dominant process by which ^Be nuclei decay in the solar interior p3[ . The following section. 



Sec. 1^, presents a calculation of the line profile when bound electrons are captured by ^Be 



ions in the solar interior. 

The partial transition probability to undergo electron capture can be written schemati- 
cally, for a given relative flux density of electrons and ^Be ions, as 

dX = flux X dcap, (23) 
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where cxcap is the appropriate weak interaction capture cross section. More specifically, 



the partial transition probability can be written in terms of the usual |^ charged-current 
beta-decay Hamiltonian, Hp, as 

d\ = 2T,\{f\Hp\l)\H (Ei,itial - ^final) , (24) 

where the continuum wave functions are assumed to be normalized in a large, finite volume. 
Averaging over initial states, summing over final states, and integrating over all space, one 
obtains 

fiux X d{cIosssecilon) = {2^l)-'^Y.Y.\U^H|3\i)\^5^^\vi-Vf)d^qd^P7, (25) 

/ i 

where the momenta in the final state of the neutrino and of the ^Li ion are denoted by g* 

and by p^, respectively. The beta-decay Hamiltonian can be written as 

Hp = 2-^/^G (^,7a (1 + 75) ^e) {i'nla {Cy - CaI,) ^p) , (26) 



where all the symbols have their usual meanings pll^3 . 



It is convenient to carry out the calculations in the center-of-momentum coordinate 
system. Relative to the laboratory frame in which the electron has a momentum pi and the 
^Be ion has a momentum p*j, the center-of-momentum frame has a velocity V given by: 

V = M-'{p,+f,), (27) 

where 

M = me + MCBe). (28) 

The equation describing the conservation of energy can be expressed simply in terms of the 
momentum, p, in the center-of-momentum frame, where p is given by 



p = M-^ (M(^Be)p; - meP7) • (29) 

One can rewrite Eq. (|]) in the form: 
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£,„,.,. - £».. = AAf + |1 - , - 5^ = 0, (30) 



where the reduced mass /z is 

^eM(^Be) 
^ me + M(7Be)' ^ '' 

The solution of Eq. ( pO] ) yields an expression for the neutrino energy in the center-of- 

momentum frame in terms of the momentum, p, of either of the particles, i.e., 

q = q{p). (32) 

The beta-decay matrix element in Eq. ( ]25|) can be greatly simplified by using the so- 



called "normal approximation" ||4T|-^ that gives the leading term for the matrix element 
in a rapidly convergent power series expressed in terms of the small (^ 0.01) quantities 
(J-Rnucieus and p-Rnucieus, whcre -Rnucieus is the nuclear radius of ^Li. The conditions for the 
applicability of the normal approximation in stars are given in Eq. (10) of ref. |^3[ and are 
satisfied to high accuracy in the present case. 

Making the usual non-relativistic reduction and summing over the spins of the nuclei in 
the initial and final states and over the spins of the electron and the neutrino, one finds in 
the center-of-momentum frame a simple expression for the product of the relative velocity 
of the particles times the capture cross section, namely, 

''^-^''--P- 27r[l + q{p)/Mim)y ^""""^ 

where \ipe\'^ is the electron probability density (enhanced by the Coulomb attraction) aver- 
aged over the nucleus and ^ is the familiar sum of reduced matrix elements that occurs in 



allowed beta-decays ||4l|-^, i. e., 

^ = CUl)' + Cl{ay, (34) 

where (1) and (a) refer, respectively, to Fermi and Gamow- Teller matrix elements. 

In a spherical shell at temperature T that contains A^(^Be) total ^Be ions, the rate at 
which ^Be captures electrons is 
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d Rate(T) = Ar(^Be) {n{e)(r,,pV,ei) , (35) 

where the average that is indicated in Eq. (|35|) is over the thermal distributions of the 
electrons and the ^Be ions. Writing out the thermal average explicitly, one finds 

^ ' 27r \27ikTj \ 2nkT J 

X f f d^p^d^prexp [-pl/2mekT) exp (-p?/2M(^Be)fcT) 

IW) \ q (p) 
^[1 + g/MrLi)]' ^^^^ 

where (il^ is the volume of the spherical shell and n(e) and n(^Be) are, respectively, the 
local number density of electrons and of ^Be ions. Converting the integration variables to 
the center-of-momentum quantities, p and P = MV [see Eq. (pT]) and Eq. (pQ])], one obtains 
the relatively simple looking equation, 

d Rate(T) = CiT) J J d'pd'Pexp [-p" /2iikT) exp (-P^/2MkT) ^^ ^'^|^^J^^^^ , 

(37) 

where C{T) is a constant that is independent of neutrino energy. To high accuracy in 
the present application, one can approximate the electron probability density by the non- 
relativistic expression pl| , |43|l . 



2 '^Tiav -^ , ^ 

We\ = 7-^ TT , Ve = fXVrcl nie, (38) 

1 — exp [Snav^^) 

where a is the fine structure constant and Vg is the velocity of the electron in the center-of- 
momentum frame. 

The neutrino energy is measured in the laboratory frame, not in the center-of-momentum 
frame. Therefore, the rate of production of neutrinos of definite observed energies, gobs? niust 
be computed. The neutrinos that are observed experience a Doppler shift because of the 
motion, V, of the center-of-momentum frame relative to the laboratory frame. Let the z 
axis be oriented along the direction between the terrestrial detector and the core of the sun. 
Then the non-relativistic expression for the Doppler shift is 
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qohs = qc.mXp) {l - Vz,c.m.) , (39) 

which corresponds to a center-of-momentuni velocity in the z direction of 

T r _ f_z^ _ Qc.mXP) ~ Qohs (Ar\\ 

Vz,c.m. ~ -nj ~ / \ ■ V^^) 

M qc.mXP) 

It is also convenient to introduce the relative energy in the center-of-niomentum frame, i.e., 

The observed energy of the neutrino depends upon the direction in which the decaying 
nucleus is moving [via Eq. (|39D], but only on the magnitude of the relative energy (not on 
the direction of p), i.e., 

QcmXP) ~ '?cont,star +P /2^ = (Jcont.star + E. (42) 

Carrying out the integrations over the unimportant directions and converting from an inte- 
gration over relative momentum to an integration over relative energy, one obtains 

rf Rate (T)= C"(T) / dE dP,exp{-E/T) exp{-P^/2MkT) 

Jo J-oo 

xql^{E) [l + q/M ('Li)]"' [1 - exp (Snamjp)]-' . (43) 



The most important terms in Eq. (|43|) have direct physical interpretations. The function 
exp{—E/T) represents the Maxwell-Boltzman distribution of relative internal energies, E, in 
the center-of-momentum frame and is dominated by the electron kinetic energy. The function 
exp(— P|/2MA;T) describes the Maxwell-Boltzman distribution of the center-of-momentum 
velocities and is responsible for the Doppler shifts (via Eq. ^9]). Since the ^Be nucleus is 
much heavier than an electron, the center-of-momentum frame approximately coincides with 
the rest frame of the capturing ^Be ion. Thus the velocities of the ^Be nuclei essentially 
determine the Doppler shifts. The last three terms in Eq. (^) represent, respectively, the 
neutrino phase space (g^), a small correction resulting from conservation of energy that is 
associated with the recoil energy of the ^Li nucleus, and a small correction to the basically 
1/f -dependence [see Eq. (|38|) ] of the probability density of the electron near the nucleus. 



The exponent that describes the Doppler shift is 

P| M(gobs-gc.„^.(^))' 



2MkT 2kTql^{E) 

which can be written in a convenient numerical form as 

Pf 51.023 (gobs -gcont,star-^)' 



(44) 



(45) 



2MkT n (1 + i5;/g,o,t,star)' ' 

where Tg is the local temperature in the sun measured in units of 10^ K and the neutrino 
energy and the internal kinetic energy are expressed in keV. The large numerical coefficient 
in the exponent that is shown in Eq. (]45| ) forces gobs to be equal to gcont,star + E to within 
a fraction of a keV. Expressed in terms of the observed neutrino energy, Eq. (^31) can be 
rewritten as: 

d Rate (T, gobs) ^,.^^ /"^ .^ / ^ ,^^ / -51.023 (gobs - gcont,star - eY\ 



roo 

C'{T) / dE exp(-E/T)exp 
Jo 



dqohs ^ ' Jo ^' ' ' '' \ Tq (1 + -E/gcont,star)^ / 

xql^{E) [l + q/M (^Li)] ' [1 - exp {-87famjp)r' , (46) 

which is the principal result of this section. 

The normalized spectrum of neutrino energies due to electron capture from continuum 
orbits within a spherical shell at temperature T is 

Q , frp ^ d Rate (T, gobs) /dq^^s ,.„. 

Spectrum^ont [T, gobs = ^^ , ^7-^ . (47) 

The average over spherical shells at different temperatures will be described in Sec. 0, where 
use will be made of detailed solar models. 

IV. CAPTURE FROM BOUND ELECTRON ORBITS 

A small, but significant fraction of electron captures in the sun occur from bound orbits. 
This section describes the calculation of the neutrino energy profile for bound electron cap- 
tures. The rate of bound-capture in the sun was first evaluated by Iben and his collaborators 



2JI in an elegant paper which carried out a variational-principle calculation of the binding 
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energies and eigenfunctions in the presence of the solar plasma and which also presented 
formulae for the fractional occupation of different bound atomic levels as a function of the 
ambient physical variables. 

The normalized energy profile resulting from bound capture has the simple form 

.2 \ 1/2 



q f . ^^_ -1 (MCBeJc 

DpeCT;rum{3Q^j^^ Wobs, -L ) — Q'bound.star o„LT^ 



Zttki 



M(^Be)c^ /gobs -gbound,star\ 



-^Kl y Abound, star / 

(48) 



where gbound,star is given in Eq. (pO]) and Eq. (p2D. The spectrum given by Eq. (|i8D is obtained 
directly from the Maxwell-Boltzman distribution for the ^Be ions by substituting for Pz the 
expression given in Eq. (^D|) for the ion momentum. For the capture of bound electrons, the 
center-of-momentum frame coincides with the rest frame of the ^Be ion, which is why the 
kinematic complications that are present for the case of continuum capture (in which both 
electrons and ^Be ions have non-zero translational velocities) are absent for bound capture. 
The fraction, /bound/(l-0 + /bound), of electron captures that occur from bound K-shell 



orbits of Be at a fixed temperature is given by the following expression ||24 



/bound(T) = (5.07/T6) Snexp {2.515aR/n) , (49) 

where Sr and a^ are quantities that result from the variational principle calculation. Bahcall 
and Moeller |35| have given a convenient formula for aji [see Eq. (0) for the definition of 



an]. They give 

aji ^ -0.431 + 2.091r - 1.481r2 + 0.401r^ , (50a) 

in terms of the dimensionless Debye-Huckel screening length, r (the Debye-Huckel screening 
length divided by the Bohr radius), where: 

r = 0.298[64r6/p(3 + X)]^/^ , (50b) 

with p and X being, respectively, the local density (in gm cm~'^) and the hydrogen mass 



fraction. The quantity Sr can be calculated from formulae given in [^ 
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1. e. 



^R = C2D-i[l + 0.435Lflexp(-0.735aR/T6)] , (50c) 

D=[i + Lr + 0.251% exp (-O.rSSa/j/Tg)] , (50d) 

Lr = 0.246 (p/x-iTe"'/') exp (2.515a^/T6) , (50e) 



where fie is the electron mean molecular weight, and |35 



Cl ^ -0.6064 + 4.859r - 5.283r^ + 1.907r=^ . (50f) 

The bound enhancement, /bound, averaged over the conditions of a standard solar model 
yields a value of /bound ~ 0.21 [^ (see also column 5 of Table 3 of the present paper). 

V. CHARACTERISTICS OF THE ENERGY PROFILE 

This section describes the characteristics of the energy profiles of the two ^Be neutrino 
lines that are shown in Fig. 1. The calculated profiles have been averaged over the physical 
parameters of detailed solar models. The discussion makes use of the energy spectra that 
were computed for fixed temperatures in Sec. ^ for electron captures from continuum 



orbits [see especially Eq. (^31) and Eq. (Wn) ] and in Sec. |3 for captures from bound orbits 
[see especially Eq. (|8]) and Eq. (||)]. 

The energy spectra are given numerically in Tables I and II. The most striking charac- 
teristic of the line profiles is their asymmetry. Figures 2 and 3 display the calculated line 
profiles. The characteristic shapes, above and below the energy at which the probability 
for neutrino emission is a maximum, will be explained physically and derived analytically 
(approximately) in the following section. Sec. |VJ 

The basic properties of the line profile that are computed numerically in this section are 
the shift in average neutrino energy. A, relative to the laboratory energy of the neutrino 
line, the full-width-at-half-maximum (FWHM) of the line profile, the half-width at half- 
maximum for energies above {W^) and below (W^) the peak of the line profile, and the 
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first, second, and third moments of the energy distributions. For all the solar models that 
are described in Table III, these characteristics of the energy profile are given in Table V 
for both ground-state capture (862 keV, see Fig. 1) and excited-state capture (384 keV, see 
Fig. 1). 

A. Solar Models 

The principal characteristics of the four solar models |1I^JT5| , |4^ that were used to compute 
the averaged line profiles are summarized in Table III. The preferred solar model, which is 
listed first in both Table III and Table V, was computed by Bahcall and Pinsonneault 



| 15| and is the only detailed solar interior model published to date that includes helium 
diffusion. For comparison, I calculate the line profile with a model labeled (No diffusion) 
that was computed |T5[ with the same input data that was used in calculating the preferred 
model, except that the "No diffusion" model does not include helium diffusion. As a test 
of the robustness of the results, I have also calculated the energy profiles with the aid 
of standard solar models computed earlier, in 1988 ||n[ and in 1982 |^^, which used less 
accurate input data. Note that the number of spherical shells used in the interior region in 
which ^Be neutrinos are produced increases monotonically from 19 shells used in the 1982 
solar model to 127 shells used in the 1992 models. The results given below show that there 
are no significant differences among the line profiles computed with the three high-resolution 
models (with interior shells > 87). 

The solar models listed in Table III span a decade of state-of-the-art solar research. All 
four models have the same central temperature to within ±0.5%, 

Tc= (15.58 ±0.08) X 10^ K. (51) 

The central temperature for the three more-precise models varies over a range of ±0.3%. The 
less accurate 1982 model has a central temperature that differs by 0.8% from the average 
modern value. 
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The average temperature, < T >7Be, computed by weighting the temperature in each 
spherical shell according to the ^Be neutrino flux produced in that shell, 

<T>7Be=E^<^CBe,T)T, (52) 

T 

is the same to within ±0.3% for the three solar models from 1988 and 1992, i.e., 

(T)7Be = (14.10 ± 0.04) X 10^ K. (53) 

The 1982 model has a ^Be-weighted temperature of (T)?^^ = 13.6 x 10^ K, which differs by 
about 3.5% from the more accurate later models, presumably because of the small number 
of shells (1^) (used in the earlier calculation. The variation over the past decade in the 
computed total ^Be neutrino flux is ±6% ; the decadal variation in the computed ^B neutrino 
flux is also ±6%. 

The relative number of electron captures that occur from bound orbits, /bound, is robustly 
determined by the solar models. From Table III, 

/bound = 0.221 ±0.006. (54) 



This result is in good agreement with the 1969 value |]3S| of /bound = 0.21. 

Table IV lists the central temperatures for ten other recently-published solar models that 
were calculated by different groups using different computer codes and different input pa- 



rameters [^], as well as the 1988 and 1992 (no-diffusion) solar models described in Table III. 
None of the solar models listed in Table IV includes helium diffusion. The heterogeneous 
set of models referred to in Table IV were derived for a variety of different applications, 
most of which were not directly related to solar neutrinos. The applications often did not 
require the highest-obtainable accuracy for the solar interior conditions. The precision with 
which the models were constructed and the accuracy of the input data varies from model 
to model. In many cases, the authors did not use the best-available radiative opacities and 
nuclear reaction rates. The central temperatures given by the variety of solar models listed 
in Table IV can be summarized by the relation 
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Tcentrai(no diffusioii) = 15.55(1 ± 0.01) X 10^ K. (55) 

All of the central temperatures in Table IV would be increased by ^ 0.1 x 10^ K if diffusion 
were included. 

Taken together, the results shown in Table III and Table IV demonstrate that the central 
temperature of the sun is determined to ±1% even without requiring unusual precision and 
accuracy in the calculations. In what follows, I shall only make use of the four solar models 
listed in Table III, since these are the only models, with which I am acquainted, whose 
characteristics are published in sufficient detail to permit precise calculations of the ^Be line 
proffie. 

B. Properties of the Energy Spectrum 

The neutrino energy spectrum that is produced in a spherical shell of the sun that is 
at a specified temperature, T, is obtained by adding the normalized energy spectra for 
the capture of electrons from continuum orbits (computed in Sec. |IIip and the normalized 
spectrum from bound orbits (computed in Sec. P^ ). The relative contributions are weighted 
by the factor /bound that was defined in Eq. (^9]). Thus 



Q , ^rp X _ Spectrum^^^t (T, gobs) + /bound (r)Spectrumb„„^d JT, gobs) , . 

Dpecxrum [i , gobsj — i ; 7 7777^ • [oO) 

-L + /bound I -f ) 

The neutrino spectrum for the entire sun predicted by a particular solar model is the weighted 
average of Spectrum (T, gobs), weighted with respect to the ^Be flux that is produced at each 
temperature. Therefore, the neutrino spectrum predicted by a given solar model is 



Spectrum,oi^^ (gobs) = J2^^ {^^^ ^) Spectrum (T, gobs) , (57) 

T 

where the ^Be neutrino flux is normalized to unity when integrated over the whole star, 

5^(i(/)(^Be,T) = 1.000. (58) 
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The values of the weighting factors, dcf) (^Be, T), are given in pubhshed tables |T^JT5|J4^ for 



all four of the models used here. The numerical values for SpectrumgQij^j.(gobs) are given in 
Tables I and II. 

The predicted probability distribution for the neutrino energy spectrum has a maximum 
at a well-defined energy, gpeak, whose location is obvious in Fig. 2 and in Fig. 3. The peak 
of the probability distribution exceeds by a small amount, 5q [see Eq. (0) and Eq. (|T^)], 
the neutrino energy, gcont.star, that corresponds to capturing an electron from a continuum 
orbit with zero kinetic energy , i. e., 

^Q — 9peak ■" ?cont,star- (59) 

For the most accurate available standard solar model, which is labeled "Bahcall-Pinsonneault 
1992 (Helium diffusion)" in Table III and Table V, the shift 5q for the (more energetic) 
ground-state transition is 

5q (g.s.) = 0.23 keV, (60) 

and is 

5q (ex.s.) = 0.04 keV (61) 

for the (less energetic) excited-state transition. Essentially identical values of 5q are obtained 
with the other solar models used here (see the last column of Table V). The value of gpeak 
can be calculated using the last column of Table V and the relation given in Eq. (0), i.e., 

Q'peak = Q'cont.star + ^Q- (62) 

The value of q-peak for ground-state transitions is 

gpeak (g.s.) = (862.27 ± 0.01) keV , (63) 

0.43 keV larger than the laboratory decay energy (see Sec. |11 Aj). For excited-state decays. 



the energy peak occurs at 
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gpeak (g.s.) = (384.47 ± 0.01) keV, (64) 

0.19 keV larger than the laboratory decay energy (Sec. [11 A| ). 

It is convenient to define the n-th moment of the solar energy spectrum by the relation 

roo 

{{q - gpeak)") = dq Spectrum,„i^,(g) {q - gpeak)" • (65) 

The moments are computed about the energy gpeak at which the probability distribution 
peaks. 

The shift in average neutrino energy, A, from the laboratory value to the solar interior 
value is 



A = i^q- giab), (66) 

or, using the definition of the spectrum moments given in Eq. (|65D, 

A = (g — gpeak) + (O'peak ^ ?cont,star) + (^cont.star ^ Q'Lab) • (67) 



The first two terms in Eq. (^7p are given in columns 3 and 9 of Table V. The last term 
can be computed from Eq. (|), Eq. (^, Eq. (|T2|) and Eq. (^Tj) of Sec. y. The value of 
(O'cont.star — fjLab) is 0.20 kcV for grouud-statc transition and 0.15 keV for excited-state tran- 
sition. Combining the results given in Table V, the shift for the ground-state transition is, 
for all three of the modern (1988-1992) models, 

A (g.s.) = 1.29 keV, (68) 

with a spread of only ±0.5%. The earlier (1982) model yields A(g.s.) = 1.23 keV. The shift 
for the excited-state transition is 

A (ex.s.) = 1.24 keV (69) 

for the three modern models with a spread of only ±0.5%. The earlier (1982) model gives 
a 4% smaller value, 1.18 keV. 
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Table V also presents some of the other calculated characteristics of the energy profile. 
The full-width-at-half maximum of the profile is denoted by FWHM and is listed in column 
4 of Table V; the half-width on the low-energy side of the peak (column 5) is denoted by 
W- and the half-width on the high-energy side (column 6) is denoted by W^+. The most 
accurate values for the ground-state transition are 

FWHM (g.s.) = 1.63 keV ; W^-(g.s.) = 0.56 keV ; l^+(g.s.) = 1.07 keV. (70) 

The corresponding values for the excited-state transition are 

FWHM (ex.s.) = 0.97 keV ; l^_(ex.s.) = 0.24 keV ; l^+(ex.s.) = 0.73 keV . (71) 

Figures 2 and 3 are remarkably asymmetric. The degree of asymmetry can be quantified 
by taking the appropriate dimensionless ratio of the third and the second moments that is 
known as the skewness. For our case, the skewness of the profile is given by 

Skewness = ^^^ ~ ^p'^""^/ . (72) 

The values of the skewness that are computed from columns 7 and 8 of Table V are, respec- 
tively, 

Skewness (g.s.) = 5.7 ; Skewness (ex.s.) = 5.5 . (73) 

The skewness of any symmetric distribution, such as the normal distribution, is equal to 
zero. 

The shape of the energy profile can be well-described by simple analytic functions. For 
neutrino energies less than the peak energy, the profile is essentially Gaussian, i. e., 

Spectrum^^i^^ (gobs) = iVexp [- (gobs - ^peak)^ /2w^] , gobs < ^peak- (74) 

The high-energy tail is well described by a Boltzmann distribution, i. e., 

Spectrum^oiar (?obs) = A^ exp [- (gobs - gpcak) /Teff] , gobs > gpeak- (75) 
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In Eq. ( [f4D and Eq. ([fSD, the quantity A^ is a normalization factor. For ground-state capture, 
the effective width, w, to be used in Eq. (|7iD is approximately 

w{g.s.) = 0.48 keV . (76a) 

The corresponding value for excited-state capture is 

w(ex.s.) = 0.20 keV . (76b) 

Both transitions are well-described by a single effective temperature, 

fcTeff = 1.31 ± 0.02 keV , (77a) 



Teff ^ 15.1 X 10*^ K. (77b) 

The uncertainty indicated in the value of the effective temperature, T^s, reflects the fact 
that Eq. (W^ is only approximate and the best-fitting value of Teff varies somewhat with 



neutrino energy. 

VI. APPROXIMATE ANALYTIC DERIVATIONS OF THE ENERGY SHIFT 

AND ENERGY PROFILE 

This section gives approximate analytic derivations of the principal characteristics of 
the neutrino energy spectrum from ^Be electron capture in the sun. The purpose of this 
discussion is to provide insight regarding the physical origins of the profiles that are shown 
in Fig. 2 and Fig. 3. For simplicity, I will usually concentrate on the profile created at a 
fixed temperature, Tg, ignoring the average over the temperature distribution of the center 
of the sun. In some cases, I will also ignore the small (20%) contribution to the energy 
profile that arises from captures from bound states (cf. discussion in Sec. ^V]). In all cases, 
I will omit small terms of order g/M(^Li), the ratio of the neutrino energy, q, to the mass 
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of the ^Li nucleus, and terms of order, i?/g, the ratio of the relative kinetic energy,£', of the 
electron and of the nucleus to the neutrino energy. 

These approximations permit the isolation of the principal physical processes that deter- 
mine the neutrino energy profile and allow analytic calculations to be made that reproduce 
the general features of the detailed numerical results. 

I begin by deriving approximate expressions for the energy half- widths, W- and W+, 
and then obtain an analytic expression for the energy shift, A. 

A. Energy Half- Widths 

The low-energy side of the profile, gobs < Q'peak, is produced by ^Be nuclei that are moving 
away from the observer and which capture electrons with small or zero kinetic energies (i. e., 
low- velocity electrons from the continuum or electrons in bound orbits). For such captures, 
the reaction rate is given approximately by [see Eq. (^3|), Eq. (^51) , and Eq. (^Sf )] 

51 



d Rate ~ const x exp 






(78) 



Replacing in Eq. (|78D the expression (gobs — gcont.star — ^g) by W-, the half-width-at-half 
maximum on the low-energy side, one has 

W-^{ln2{n)/5lf\ (79) 

Evaluating Eq. (^) for a characteristic ^Be temperature of 14 x lO^K (see discussion of solar 
models in Sec. [VA|) yields 

W- ^ 0.44 keV. (80) 

This result is in satisfactory agreement with the low-energy width that is calculated numer- 
ically, W^ = 0.56 keV (see Table V). 

The high-energy side of the profile, gobs > gpeak, is produced by electrons and ^Be nuclei 
with significant internal kinetic energies, i? [see Eq. (^3|)]. Since the product of the separate 
Boltzmann distributions for the electrons and for the ^Be ions is a Boltzmann distribution 
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in the center-of-inomentum frame, the high-energy part of the profile has the approximate 
shape 

d Rate ~ const x exp{-E/T) , (81) 

where 

E = pV2fi, (82) 



and the relative momentum, p, is defined by Eq. (pQl). One can define an effective tempera- 
ture, Tofj, for the high-energy tail of the spectrum by fitting at two energies the expression 
given in Eq. (^) to the numerically-computed energy spectra that are given in Table I and 
Table II (and shown in Fig.3 and Fig.4). Symbolically, 

gobs (2) -gobs(l) 



^eflf 



(83) 



In [Spectrumobs(gobs(l))/Spectrumobs(gobs(2))] ' 
where gobs(l)5 Q'obs(2) are typical energies in the hne profile. The value of Tefj depends slightly 
on the choice of qohsi^), Q'obs(2). For both transitions shown in Fig. 1, 

Teff ^ 1.31 keV ^ 15 X 10^ K, (84) 

which translates into a half-width-at-half-maximum of 

W+ ^ (ln2)Teff ^ 0.91 keV. (85) 

The value for W^ given in Eq. (^5]) is approximately halfway between the accurate values 
given in Table V for the ground-state decay (1.07 keV) and the excited-state decay (0.73 
keV). 

The reason that the effective temperature given in Eq. ( |8^ is slightly larger than the value 
obtained from the solar model and given in Eq. ( |5^ ) is that Eq. (p3| ) is only approximate. 
A more accurate (but less transparent) determination of T^s can be obtained by taking the 
ratio of integrand in Eq. (^) at two neutrino energies and setting that ratio equal to the ratio 
of the spectrum probabilities. This more accurate relation should be used in interpreting 
future experiments. 
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B. The Energy Shift: A 

The energy shift, A, can be written as the sum of two terms, (gobs — Q'cont,star) plus 
(gcont.star " ^lab) [cf- Eq. (§7^]. FoT the ground-state transition, (gcont.star - giab) is 0.20 keV 
and is 0.15 keV for the excited state transition [see discussion in Sec. ^ and in Sec. |V B| 
following Eq. (|57D]. The dominant term in the energy shift, (gobs — ^cont.star), is equal (to 
numerical accuracy, which is better than 1%) to the same value, 1.09 keV, for both ground- 
state and excited-state transitions. This numerical equality can be established as follows. 
The quantity (gobs - gcont.star) is equal to the sum of two terms, (gobs - "Jpeak) + Sq, both 
of which are given in Table V. Adding the two terms in Table V gives, for all of the solar 
models, essentially identical values for (gobs — gcont,star)- 

At first sight, it is surprising that the energy shift relative to gcont.star, (^obs - gcont,star), 
is independent of the value of gcont,star- However, a simple analytic argument, given below, 
shows that, at a specified temperature T, the shift (gobs — Q'cont,star)T is, to an excellent 
approximation, the mean kinetic energy of the electron and the ^Be nucleus that interact to 
produce the electron capture reaction. The mean kinetic energy is, of course, independent 

OI gcontjStar- 

The average value of the shift relative to gcont.star is 

(^obs ~ Q'cont,star)T = / '^'?obs ((Jobs " "Jcont.star) Spectrum^(gobs), (86) 

Jo 
which can be rewritten approximately using Eq. (^31) and Eq. ( ^4|) as 

_ f fJTT Y^^ f^ 

(gobs - gcont,star)T = T M — / dE exp{-E/T) 

\gcont,star/ "'O 

X / dqohs (gobs - gcont.star) exp -/? (gobs " gcont,star " E) /gcont,star • (§7) 



The key fact that simplifies Eq. (|87D and results in the equality of the energies for ground- 
state and excited-state transitions is that the quantity, /3, which appears in the second 
exponent in Eq. (|57|), is very large. Eq. (^) and Eq. (^) show that P is just the ratio of 
the nuclear rest mass energy (~ 7 GeV) to the ambient temperature (~ 1 keV), i. e., 
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.^i^»i. (SS) 



In the limit of large (3, 



(/5/7r)^/2 g^p (-px^) — . 6{x), (89) 

which greatly simplifies Eq. (^). The effect of the delta function is to enforce energy 
conservation independent of the Doppler shifts, i. e., 

Q'obs — Q'cont,star = E. (90) 

Carrying out the integration over the delta function gives 

POO 

(gob. - gcont,star)r = T'^ / dEE cxp {-E/T) = T (91) 

Jo 

Therefore, 

(gobs - gcont,star)0 = J dTT d<P [' Be , t) = 1.2 keV, (92) 

which is, as promised, independent of gcont,star- 

VII. OTHER EFFECTS 

This section shows that the effects on the neutrino energy profile of electrostatic screening 
energy (Sec. [VII A|) , of the gravitational redshift (Sec. [VII B| ), and of collisional broadening 
(Sec. [VII C] ) are much smaller than the dominant terms (calculated in Sec. fT^Sec. |V|) that 
arise directly from the Maxwellian energy distributions of the electrons and of the ions. 

A. Electrostatic Screening Energy 

There is a difference in electrostatic screening energy between the initial and final states 
of ^Be continuum electron capture. In the initial state, a ^Be nucleus (charge z = 4e) and a 
nearby electron jointly interact electrostatically with a surrounding charge cloud of electrons 
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and ions that screen the nucleus. After the electron is captured, only the nuclear Coulomb 
potential corresponding to a charge z = 3e interacts with the screening cloud. The difference 
in electrostatic screening energy is contributed partly to the energy of the emitted neutrino. 
The electrostatic screening energy, 6Ecoui, can be estimated by generalizing the familiar 
Debye-Huckel calculation of the charge distribution and the potential associated with the 
monopole field of a single ion. The more general case consists of a monopole field (from 
the capturing ion) plus a dipole field (from the ion plus the electron that is about to be 
captured). Solving Poisson's equation for the monopole plus dipole fields, one obtains the 
potential, ^before (j), before the electron capture occurs. 



Wore [r) = Znet \ exp(-/tr) 



1 + (/tr) 



-1 



cos 6, 



(93) 



where k is the inverse of the Debye-Huckel screening length: 



v^ = Ane^T-^^^zfri,. 



(94) 



The first term on the right hand side of Eq. ( PHj ) is the standard Debye-Huckel potential for 
a point charge, ^net, surrounded by a screening cloud, where the net (point) charge in the 
present case is +3e. The second term is the dipole solution in which p is the dipole moment 
produced by the ^Be nucleus and the (to-be-captured) electron. In the usual approxima- 
tion that the screening potential is smaller than the thermal energy, the charge density 
surrounding the ^Be nucleus and the electron is 



K 



Pbefore (t) = -"— ^before (^)- 



(95) 



The Coulomb energy both before and after the electron capture can be evaluated from 
the relation 



^coui = - / d-'rpirmr). 



(96) 



In evaluating the initial Coulomb energy, the potential and the charge distribution are taken 
from Eq. (^) and Eq. (0) . One may approximate the final Coulomb energy by assuming 
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that the density distribution of the charge cloud is unchanged immediately after the electron 
capture occurs, but that the final potential only includes a pure monopole term (3 e/r) due 
to the final nucleus plus the potential due to the screening charges. The potential due to 
the screening charges can be taken from the following defining relation 

iP T) 

0beforc (r) = ^screening (^) H h —COS 6. (97) 

The difference in Coulomb energies, SEconi, obtained in this way is 

where Xd is the De-Broglie wave length of the electron and Rbu = k~^ is the Debye-Huckel 
screening length. In obtaining Eq. (|98|), I evaluated the dipole moment by assuming that 
the electron and the ^Be nucleus are separated by a distance Ad with the origin at the center 
of charge. This yields p = ScAd/S. Inserting typical solar-interior values into Eq. (|98D, one 
finds 

^^approx. Coul ~ "0.004 keV. (99) 

One can obtain a conservative upper limit to the change in electrostatic screening energy 
by assuming that the electron is far away from the nucleus when the electron is captured. 
In this extreme limit, one can approximate the change in Coulomb energy by 

^coui = ^J d^rp{r)z=4 [0(r )z=4 - 0(r )z=3] , (100) 



where in writing Eq. (|100|) it was also assumed that the particles in the screening cloud 



around the ^Be nucleus do not have time to move during the electron capture. These 
assumptions yield 

.^^extreme Coul ^ - -— ^ "0.065 keV. (101) 



Some fraction, e < 1, of the Coulomb energy difference represented by Eq. ( PP| ) or Eq. 
( |101|) should be added to the neutrino energies calculated in the previous sections. This 
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addition is estimated to lie between 0.004 e keV and 0.065e keV, which is a rather small 
effect. However, a quantum mechanical calculation of the change in electrostatic screening 
energy would be of interest. 

B. Gravitational Redshift 

Each neutrino energy is shifted by an amount Ag, 

A,(G.R.) ^ Z^^^,. (102) 

due to the general relativistic redshift. In the region in which the neutrinos are produced, 



the mean shift Ag/g ^ 10 ^ ||28|. Therefore, the average change in the neutrino energy is 
only 

Ag(G.R.) ^ -0.009 keV. (103) 

This energy redshift is two orders of magnitude smaller than the dominant processes that 
are calculated in Sees. [m]^VI|. The dispersion in the gravitational energy shift, which would 
contribute to the broadening of the line, is an order of magnitude smaller, ~ 0.001 keV. 

C. CoUisional broadening 

The effect of collisions on the coherence length for neutrino oscillations was first discussed 



by Nussinov |^ in a beautifully-original paper(that also quantified the small probabilities for 



obtaining a large reduction in the electron-neutrino flux at earth due to vacuum oscillations). 



Loeb [0 was the first to carry out a detailed calculation of the collisional broadening of a 
solar neutrino line (see also |^8|]). In the discussion below, I follow the treatment of Loeb. 

Let -P(g) be the probability that a neutrino is emitted with the energy q in the presence 
of the solar plasma and let go be the energy that would be emitted if the collisional frequency 
were set equal to zero. The probability distribution P{q) has the usual Lorentz shape. 
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m = , ^"'t^ , , (104) 

(9-90) +qr 

where 

q, = {2tit)-\ (105) 



and r is the coherence time for neutrino emission |^. The coherence time denotes the 
period over which the emitting system looses its phase due to colhsions with the background 
plasma particles. 

Loeb 1^ shows that the collision time for ^Be ions in the solar interior is about lO^^s^^ 



and that the coherence time, r, is 

r?a5x 10"^^ sec. (106) 

The width, g,-, is therefore very small, 

Qr ^ 0.013 keV, (107) 

and may be neglected in the present context. 

VIII. NEW PHYSICS 

The two most popular mechanisms for explaining the solar neutrino problem via new 
physics are vacuum neutrino oscillations, first discussed in this connection by Gribov and 
Pontecorvo |l^ in an elegant and epochal paper, and matter-enhanced neutrino oscillations. 



the MSW effect, a beautiful idea discovered by Wolfenstein [17] and by Mikheyev and 



Smirnov ||T8[. In Sec. |V111 A|, I present calculations of the effect of vacuum oscillations on 



the energy shift. A, of the ^Be neutrino line (solar versus laboratory) and in Sec. |V111B| , I 
investigate the effect of matter-enhanced oscillations. Finally, in Sec. |V111C| , I discuss briefly 



the effect of other suggested new physics processes on the energy shift. These discussions 
of potential (new-physics) solutions show that the physical processes considered would not 
be expected to change significantly the energy shift that is calculated by considering only 
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thermal effects in the core of the sun. The effect of the width of the ^Be hne on neutrino 



oscillations was discussed previously by Pakvasa and Pantaleone |26 



A. Vacuum Oscillations 

If vacuum oscillations occur, the (energy-dependent) probability that an electron-type 
neutrino, Ve-, that is created in the solar interior with an energy gobs in the laboratory 
frame survives as an electron-type neutrino until it reaches a terrestrial detector modifies 
the original solar neutrino energy spectrum. Thus 

Spectrum^^ (gobs) = Spectrum^oiar (?obs) |(z/e, detect |i/e, emit) |g^^^. (108) 

If the electron-type neutrino is primarily coupled to only one other neutrino type, e.g., a 
muon-type neutrino, z/^, then the energy spectrum of the daughter neutrino is 

Spectrum^^ (gobs) = Spectrum^oiar (?obs) 1 - |(z/e, detect |i/e, emit) 1^^^^^ . (109) 

The probability for an electron-type neutrino to change its flavor to a muon-type neutrino 
can be written as |]16|,E9lE0 



|(z/^,detect|z/e,emit)|q^j^^ = 1 - |(z/e, detect|i/e,emit)|g^i^^ = sin^ 29v' sin^ (gobs) , (HO) 

where 9y is the vacuum oscillation angle and 0(gobs) depends upon the neutrino energy, 
upon the differences in the square of the masses of z/^ and i/g, Am^, the distance, D, be- 
tween the terrestrial detector and the location in the sun where the neutrinos are produced. 
Numerically, 

, , n /0.86 MeV\ /Am2\ / R \ ^ ^ 

^(,„,, . ,15 X 10" [——) -^ (^ . (HI) 



As pointed out by Bahcall and Frautschi |H9[ , the survival probability must be averaged 



over the neutrino energy profile in order to calculate the effect of vacuum oscillations on the 
observed neutrino event rates. Fine-tuning of the vacuum oscillation parameters is required 
to produce large effects on the observed rates. If vacuum oscillations reduce the flux of ^B 
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electron- type neutrinos (which have a broad energy profile with a total width ~ 10 MeV) by 
a large factor (as required to explain the difference between the standard model predictions 
and the observations with the chlorine and Kamiokande detectors), then the mass difference 
squared, 5rn?, must cause 0(gobs) (evaluated at ^ 8 MeV) to be a small integer multiple of 
7r/2 [^. Thus, if vacuum oscillations are to explain the solar neutrino problem, 

Am^ (big effect;^ b) ~ IQ-^"-^ eV^. (112) 

There have recently been several careful studies |^ of the constraints on the vacuum 



oscillation parameters that are implied by the existing four solar neutrino experiments @-0 . 



Krastev and Petcov [^ summarize these results as follows: 

5 X 10"" eV^ < 5m'^ < 11.1 x 10"" eV^ (113) 

0.75 < sin2 2ey < 1.0. (114) 

Only certain combinations of 5rn? and sin^ 2By are allowed, but for convenience and to be 
conservative, I have explored the entire range in the rectangular space defined by Eq. ( |113[ ) 
and Eq. ( |114D . The allowed range of (gobs = %eak) is 



3.5 X 7r/2 < 0(862.27 keV) < 8.0 x 7r/2. (115) 

The survival probability for an electron-type neutrino does not change much over the 
1 kev width of the ^Be line. In order for the phase-angle, 0(gobs) (defined by Eq. (|111|) 
to change by an appreciable fraction of a radian, the observed energy must change by an 
amount, Agobs, that is much greater than the full width of the neutrino line. The line profile 
is not significantly affected by vacuum oscillations unless the phase-angle, 0, is chosen just 
so as to make the electron-neutrinos at the peak of the energy ^Be energy profile maximally 
mix into muon neutrinos, i. e., 0(gpeak) is chosen to be an odd integer multiple of 7r/2. 
Even in the case of maximum mixing, for which electron neutrinos are practically all fiavor- 
converted, the resulting muon neutrinos have essentially the same energy spectrum as the 
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original spectrum with which the electron neutrinos are created. Therefore, the energy 
profile that would be detected by neutrino-electron scattering is essentially the same for 
maximum mixing and for no mixing. The observed rate is decreased, of course, by mixing 
because muon neutrinos scatter off electrons about a factor of five less strongly than electron 
neutrinos. 

How can we quantify the effect of vacuum oscillations on the ^Be neutrino line profile? 
The most direct effect of oscillations on the ^Be line profile is manifested in the difference, A 
Eq. (^), between the average energy of solar-produced neutrinos versus laboratory-produced 
neutrinos. Table VI presents results of numerical calculations that have been carried out 
by averaging the effect of vacuum oscillations over the energy profile of the ^Be line that is 
computed using the Bahcall-Pinsonneault solar model with helium diffusion (see Sec. [V A| ). 
The quantity A(z/e) that is given in Table VI is 

, , / rfgobsSpectrum (gobs) (l - sin^ 2Qy sin^ (gobs)) (^obs - giab) 

A(z/e) = ^ 7 X , (116) 

/ rfgobsSpectrum (gobs) (l - sin^ 2Qy sin^ (gobs)) 

which is the energy shift that would be measured for electron- type neutrinos (e.g., in neutrino 
absorption experiments). The corresponding energy shift, A(z/^), that would be measured 
with muon-type neutrinos is 

., . _ I (igobs Spectrum (gobs) sin^ (gobs) (9obs - 9iab) .^.„s 

^ ~ /rfgobsSpectrum (gobs) sin^^ (gobs) 

For neutrino-electron scattering experiments, the energy shift that would be measured, 

(A)el.sc., is 



fJeSurvival (ue) Ag + cr^ (1 — Survival (i^e)) A^ 
'^''^'^' (JeSurvival (z^e) -|- cr^ (1 — Survival (z^e)) 



_ (^e^mvivdiyt^ej^e-r c^^yx - ^ULVivaL^iyeJJ^fM ^^gx 



I list in Table VI representative numerical results obtained for the two extreme values of 
the vacuum mixing angle, sin^ 20y = 1.00 and sin^ 2Gy = 0.75 [cf. Eq. ( |114D ]; intermediate 
choices of the vacuum mixing angle yield intermediate effects that can be approximately 
interpolated from Table VI. I do not list values when the component in question carries less 
than 10% of the flux, since these cases are not relevant for currently feasible experiments. 
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The numerical calculations show that the energy shift, A, for both electron and muon 
neutrinos is within 10% of the shift, 1.29 keV, calculated in the absence of vacuum oscillations 
for all cases in which the corresponding neutrino flux (electron or muon) is not less than 
10 % of the total flux. A change of a given sign in the energy shift for electron neutrinos, 
A(z/e), implies a change of the opposite sign in the energy shift for muon neutrinos, A(z/^). 
The calculated energy shift, < A >ei.sc., which should be used in making comparisons with 
electron-neutrino scattering experiments, is always within 3% of the no-oscillation value of 
1.29 keV. 

B. An MSW Solution 

The MSW solution, matter-induced neutrino oscillations, to the solar neutrino problem 
has been discussed by many authors |T7| , p!8| , p2[ . For matter-enhanced oscillations, the proba- 
bility of neutrino mixing within the sun depends upon neutrino energy, but not in as delicate 
a fashion as for vacuum oscillations. The dependence of the survival probability on energy 
can be represented as a smooth function of energy over the extent of the ^Be line proflle 



13| , p!7|jl8| , |22[| . The dimensionless ratio that determines the amount by which the MSW effect 
changes the line proflle is the ratio of the neutrino line width to the characteristic energy of 
the line. This ratio is very small, FWHM/2gpeak = 0.1%. 

Large suppressions of electron-type neutrinos are achieved by the MSW effect without 
flne-tuning with respect to the neutrino energy. In fact, it does not seem possible to flne-tune 
MSW solutions to an accuracy that would greatly amplify the small dimensionless ratio of 
FWHM to Qpeak and thereby affect signiflcantly the line proflle of a large fraction of the 
^Be solar neutrino flux. In the resonance condition, the neutrino energy, gres, is inversely 



proportional to the electron density,ne, at the resonance position [|13],|T^|I8| . The change in 
the electron density over the region of production of ^Be neutrinos is [|1^] Anc/ric ~ 0.25, 



while the change in neutrino energy over the line proflle is Aq/q < 0.001. If one tried to 
invent a situation in which only part of the ^Be line proflle went through resonance, the 
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fine-tuning that applied at one location would quickly be destroyed as the electron density 
changed by a tiny amount (and therefore changed the resonance condition) within the region 
in which the neutrinos are produced. 

C. Some Other New Physics Processes 

Other solutions have been proposed for the solar neutrino problem that involve new 
weak interaction physics. These other solutions include rotation of the neutrino magnetic 



moment pO|, matter-enhanced magnetic moment transitions |^, and neutrino decay ||52|| . 
The classical magnetic moment transition is independent of energy and does not affect 
the shape of the line profile. Matter-enhanced magnetic-moment transitions, like MSW 
transitions, are not fine-tuned, vary smoothly with energy, and depend upon a resonant 
electron density that varies from point to point. Therefore, the argument given above for 
the MSW effect also applies to matter-enhanced magnetic-moment transitions. Neutrino 
decay involves characteristic energies that are very large compared to the total width of the 
^Be line and is also not fine-tuned. Hence, none of these processes would change significantly 
the shape of the ^Be line profile. 

IX. THE ^LI NEUTRINO ABSORPTION CROSS SECTION 

The calculated absorption cross section for the reaction 

i/e + ^Li^^Be + e", Eth = 861.9 keV, (119) 

where v^ is produced by ^Be electron capture in the sun, depends upon the assumed energy 
profile of the solar neutrinos p7|j38[1 . Neutrinos with energies below the energy threshold, 
ii^th; for the laboratory reaction, Eq. ( |119| ), cannot be absorbed. The energy threshold of 
861.90 keV (cf. Eq. |^) falls within the line profile shown in Fig. II. The precise location of 
the threshold within the line profile determines the fraction of emitted neutrinos that can 
be absorbed by ^Li. 

41 



The average absorption cross section for solar-produced ^Be neutrinos incident on a 
laboratory detector of ^Li is 

(Spectrum^^(gobs) aabs(gobs)) ^ 19 x 10"^^ cm^ (120) 

assuming neutrinos do not change flavor after their creation. To the accuracy shown, the 



results are identical for ^Be line profiles calculated using the standard models |T5| with and 
without helium diffusion (cf. rows 1 and 2 of Table III). As usual, Eq. (|120D includes a 
correction to take account of the fact that only 89.7% of the ^Be neutrinos are produced in 
ground-state to ground-state transitions. 

The cross section given in Eq. (|120|) is almost a factor of two larger than obtained pre- 
viously p7| , |38| , which should make the contemplated radiochemical experiments somewhat 



easier than originally considered [Q. The earlier treatments neglected the difference in 
electron binding energies of solar and laboratory ^Be atoms as well as Doppler shifts of the 
^Be nuclei, and did not average over the temperature profile of the sun. In the present 
calculation, 88% of the ^Be neutrinos are above threshold for the reaction Eq. ( p.l9| ). For 
the earlier calculation |^, only about 50% were above threshold. The difference between 



a cumulative probability of 50% and a cumulative probability of 88% corresponds to an 
average energy shift of 0.85 keV for the line profile shown in Fig. 2. 

X. SUMMARY 

The temperature distribution of the solar core is expressed robustly in the neutrino 
energy profile that results from ^Be electron capture. The characteristics of the line profile — 
shown in detail in Figures 2 and 3 — are independent of uncertain details regarding solar 
models and instead reflect the global thermal properties of the solar interior. The robustness 
of the computed characteristics of the ^Be line profile derives from the well-determined 
thermal structure of the solar- model description of the interior of the sun (cf. Sec. |V A| ). 
For the most precise standard solar models computed over the past decade, 1982-1992, 
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and listed in Table III, the central temperature has varied over a total range of ±0.5%, 
Tc = (15.58 ± 0.08) X 10*^ K. The ^Be and the ^B neutrino fluxes have varied by ±6%. The 
relative number of electron captures that occur from bound orbits is 0.221 ±0.006 for all the 
solar models of Table III. For a heterogeneous set of ten recently-calculated solar models, 
listed in Table IV, that were generally not required to have the highest attainable precision, 
the central temperature varied by ±1%, Tc = (15.55 ± 0.15) x 10^ K. 

The following paragraphs summarize the principal results obtained in this paper. In 
addition. Tables II and III present numerical representations of the line profiles and Table V 
provides a concise summary of the potentially measurable characteristics that were calcu- 
lated using different solar models. The numerical values for the energy shift. A, and the 
low-energy and high-energy half-widths of the line profile, are calculated in Sec. 0. 

The shift in average neutrino energy. A, between ^Be neutrinos emitted in the sun and 
^Be neutrinos produced in the laboratory is A(g.s.) = 1.29 keV for ground-state captures and 
the physical parameters of the most accurate available solar model (which includes helium 
diffusion). Calculations with other models yield values of the energy shift for ground-state 
decays of 1.28 keV, 1.28 keV, and 1.23 keV (for a 1982 solar model with only 19 shells in the 
region in which "^Be neutrinos are produced). The shift for excited-state decays is A(ex.s.) = 
1.24 keV. These calculated energy shifts take account of the fact that approximately 88% 
of the captures involve electrons in continuum orbits and only about 22% involve electrons 
that are bound to the decaying nucleus. The atomic binding energies that are released when 
^Be nuclei capture electrons in the sun (or in the laboratory) are evaluated in Sec. ||. 

The low-energy half-width, W-, of the line profile is (0.55±0.02) keV for the ground-state 
decay (0.24 keV for the excited-state decay, see column 5 of Table V). Here W- is the half- 
width of the energy profile below the peak in the probability distribution. This low-energy 
part of the line is primarily determined by Doppler shifts caused by the thermal velocities 
of ^Be ions that are moving away from the detector located on earth. The low-energy side 
of the profile is approximately Gaussian in shape, reflecting the Doppler origin of this part 
of the energy spectrum. 
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The high-energy half- width, W^, of the hne profile is 1.07 keV for the ground-state 
decay (0.73 keV for the excited-state decay, see column 6 of Table V) and is determined 
primarily by the center-of-momentum kinetic energies of the electrons and the ^Be nuclei 
that take part in electron capture reactions. The high-energy side of the line profile is 
approximately exponential in shape, with a probability distribution that is proportional to 
exp{—qohs/Tcs) where gobs is the neutrino energy that is observed in the laboratory and 
Teff ~ 1.31 keV (15.1 x 10^ K). The exponential side of the probability distribution results 
from an average (in the solar model) over the different exponential distributions of center- 
of-momentum energies (electrons and ions) that apply at each solar radius. 

The principal characteristics of the line profile are derived by approximate analytic cal- 
culations in Sec. 0, calculations that elucidate the physical origins of the various effects. 



The most remarkable result obtained in Sec. VI is that the energy shift. A, is essentially 



identical for the ground-state and the excited-state captures (cf. Fig. 1). The average shift 
is shown in Sec. |V^ to be independent of the typical emitted neutrino energy (384 keV or 
862 keV) because the rest-mass energy of a ^Be nucleus is much larger than the solar thermal 
energies (see discussion in Sec. |V1B| ). 



The energy shift is, to a good approximation, equal to the average temperature of the 
solar interior weighted by the fraction of ^Be neutrinos that are produced at each tempera- 
ture, i. e. /q dTTd(f)(jBe, T)/ Jq dTdcpijBe, T), where d(f)(^Be, T) is the flux of '^Be neutrinos 
produced at the local temperature T [cf. Eq. (0)]. The ^Be neutrinos are produced in 
the inner few percent of the solar mass [|l^, essentially all in the region (0.04 ± 0.03) Mq. 



Therefore, a measurement of the energy shift is a measurement of the central temperature 
distribution of the sun. 

The most striking aspect of the computed energy profile is the asymmetry between the 
Gaussian low-energy side and the exponential high-energy side (cf. Figs. 2 and 3). Doppler 
shifts caused by the thermal velocities of the ^Be nuclei produce a symmetric, Gaussian 
contribution to the line broadening, which determines the shape of the energy profile at 
energies below the peak. The higher-energy part of the profile is determined by the center- 
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of-inoinentum kinetic energies. The exponential distribution of kinetic energies produces an 
exponential tail for large neutrino energies. The positive-definite character of the kinetic 
energies is responsible for the asymmetry. 

Electrostatic screening, particle collisions, and gravitational redshifts all contribute to 
the line broadening, but their effects are much smaller than the effect caused by thermal 



broadening (Sec. [VI I]) . 

Vacuum neutrino oscillations can be fine-tuned to produce maximal mixing of neutrino 
flavors near the peak of the ^Be neutrino line. But, the energy shift of the dominant neutrino 
survivor is always close to the unmixed value of A = 1.29 keV (see Table VI). The invariance 
of the line shape results from the fact that the oscillation phase changes significantly only 
over an energy range that is much larger than the line shape (see Sec. [VIII A|) . 

The energy profile of the ^Be neutrino line should be taken into account (using Tables II 
and III) in precise calculations of what is to be expected from vacuum neutrino oscillations. 
It has become standard in calculations of the effects of vacuum oscillations to take account 
of the variation of the distance between the point of creation of the neutrinos and the point 
of detection. The variation in the point of creation corresponds to a phase-change of order 
10~^, since the ratio of the solar radius to the earth-sun distance is about 0.005 and the ^Be 
neutrinos are produced in a region of about ±0.025-^0. On the other hand, the width of 
the ^Be line profile is about 10~'^ of the average ^Be neutrino energy. Therefore, the change 
in phase due to the energy width of the neutrino line is about an order of magnitude larger 
than the phase-change caused by averaging over the region of production. 

The MSW solutions that are consistent with existing solar neutrino experiments vary 
smoothly with energy and are not fine-tuned. The variation of the resonant electron density 
within the sun prevents fine-tuning of the solution over a small energy range like the width 
of the ^Be neutrino line. For MSW solutions, the fractional change in the electron-type 
neutrino survivability over the energy profile of the ^Be line is small. Other solutions of 
the solar neutrino problem that involve new physics like rotation of the neutrino magnetic 
moment, matter-enhanced magnetic-moment transitions, and neutrino decay are also not 

45 



expected to change significantly the shape of the ^Be line profile. 

The energy shift and the shape of the ^Be neutrino energy profile affect significantly the 
computed value of the absorption cross section for ^Be solar neutrinos incident on a ^Li 
detector (see Section (TXI) . 
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FIGURES 
FIG. 1. The ^Be Decay Scheme. The laboratory decay scheme for ^Be is shown. The neutrino 
energies emitted in the ground-state to ground-state decay (branching ratio: 89.7%) and in the 
ground-state to excited-state decay (branching ratio 10.3%) are denoted, respectively, by gLab(g-s.) 
and (/Lab(ex.s.). Details of the nuclear physics properties are summarized in [27|. 



FIG. 2. The Energy Profile for the 862 keV line. The probability for the emission of a neutrino 
with energy f/obs in the laboratory frame is shown as a function of (7obs — Q\ah, where the peak in the 
probability distribution is Qpcak = 862.27 keV, 0.43 keV larger than the laboratory decay energy, 
Qiab = 861.84 keV. The probability distribution was computed by averaging Eq. (H3),Eq. (|48|), Eq. 



3),Eq. (pQ), Eq. (56), and Eq. ( pTD over the Bahcall-Pinsonneault standard solar model with 
helium diffusion. 

FIG. 3. The Energy Profile for the 384 keV line. This figure was computed in the same way as 
Fig. 2. The peak of this distribution lies at Qpeak = 384.47 keV, 0.19 keV larger than the laboratory 
decay energy, giab = 384.28 keV. 
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TABLE I. The Energy Profile of 
ergy is measured relative to the ener, 

-P(g'obs - g'pcak) = Spectrumsoiar(9obs 
emitted between q ±0.1 keV. 



TABLES 

the 862.0 keV ^Be solar neutrino line. The neutrino en- 

gy of the peak, which occurs at g'pcak = 862.27 keV. Here 

is the probability that a neutrino of energy (7obs will be 



^obs 'Zpcak 


PiQohs - Qpcak) 


9obs 9pcak 


P(.qohs - gpcak) 


9obs 9pcak 


^(^obs - ^pcak) 


-2.100 


0.00000 


1.000 


0.26096 


4.100 


0.02526 


-2.000 


0.00000 


1.100 


0.24204 


4.200 


0.02341 


-L900 


0.00001 


1.200 


0.22456 


4.300 


0.02169 


-L800 


0.00004 


1.300 


0.20838 


4.400 


0.02010 


-L700 


0.00012 


1.400 


0.19339 


4.500 


0.01863 


-L600 


0.00032 


1.500 


0.17946 


4.600 


0.01726 


-L500 


0.00082 


1.600 


0.16655 


4.700 


0.01600 


-1.400 


0.00197 


1.700 


0.15455 


4.800 


0.01482 


-1.300 


0.00442 


1.800 


0.14342 


4.900 


0.01373 


-1.200 


0.00928 


1.900 


0.13308 


5.000 


0.01272 


-1.100 


0.01825 


2.000 


0.12347 


5.100 


0.01179 


-1.000 


0.03359 


2.100 


0.11455 


5.200 


0.01092 


-0.900 


0.05787 


2.200 


0.10627 


5.300 


0.01012 


-0.800 


0.09338 


2.300 


0.09858 


5.400 


0.00937 


-0.700 


0.14115 


2.400 


0.09145 


5.500 


0.00868 


-0.600 


0.20014 
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0.08482 
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0.00805 
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0.26666 
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0.33448 
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0.07295 
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0.00690 


-0.300 


0.39629 
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0.06765 
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0.00592 


-0.100 
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0.000 


0.48933 


3.100 


0.05394 
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0.00508 


0.100 
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0.05001 
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0.46807 
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0.00436 
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0.00321 
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0.32846 


3.800 


0.03173 


6.900 


0.00297 


0.800 


0.30396 


3.900 


0.02941 


7.000 


0.00275 


0.900 


0.28153 


4.000 


0.02726 
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TABLE II. The Energy Profile of the 384.5 keV '^Be solar neutrino line. The neutrino en- 
ergy is measured relative to the energy of the peak which occurs at Qpcak = 384.47 keV. Here 
Pilohs — Qpeak) = SpectrumsQig^j.(gobs) is the probability that a neutrino of energy gobs will be 
emitted between q ± 0.1 keV. 



Qohs 'Zpcak 


P{qohs - Qpcak) 


Qohs 9pcak 


P{qohs - 9pcak) 


Qohs Qpeak 


Piqohs - Qpcsk) 


-1.000 


0.00000 
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0.16166 


4.600 


0.01956 


-0.900 


0.00000 


1.900 


0.15004 


4.700 


0.01812 
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0.13924 
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0.01679 
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0.00022 
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0.00187 
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0.01442 
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0.01152 
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0.01336 
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0.10320 
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0.01238 


-0.300 
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0.35442 


2.600 


0.08881 
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0.000 
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0.00846 


0.200 
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0.400 


0.46715 
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0.05652 
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0.00673 


0.500 


0.42841 


3.300 


0.05241 


6.100 


0.00623 


0.600 


0.39618 


3.400 


0.04859 


6.200 


0.00577 


0.700 


0.36720 


3.500 


0.04505 


6.300 


0.00535 


0.800 


0.34060 


3.600 


0.04177 


6.400 


0.00495 


0.900 


0.31603 


3.700 


0.03872 


6.500 


0.00459 


1.000 


0.29330 


3.800 


0.03590 


6.600 


0.00425 


1.100 


0.27224 


3.900 


0.03328 


6.700 


0.00394 


1.200 


0.25271 


4.000 


0.03085 


6.800 


0.00365 


1.300 


0.23460 


4.100 


0.02859 


6.900 


0.00338 


1.400 


0.21778 


4.200 


0.02650 


7.000 


0.00313 


1.500 


0.20216 


4.300 


0.02457 


7.100 


0.00290 


1.600 


0.18765 


4.400 


0.02277 


7.200 


0.00269 


1.700 


0.17418 


4.500 


0.02110 
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TABLE III. Some characteristics of four standard solar models. 



Model 


(T 


Shells 
> 6 X 10'^ 


K) 


-'central 

(10** K) 


(T> 
(Weight: ^Be) 


Bound Fraction 
'Be decays) 


<^(^Bc) 
(IQS cm-2s-i) 


(10'^ 


>(»B) 

' cm ^s 1) 


Bahcall-Pinsonneault 1992 
(Helim diffusion) 




127 




15.67 




14.13 


0.215 


4.9 




5.7 


Bahcall-Pinsonneault 1992 
(No diffusion) 




127 




15.57 




14.06 


0.217 


4.6 




5.1 


Bahcall-Ulrich 1988 
(No diffusion) 




87 




15.64 




14.07 


0.219 


4.7 




5.8 


Bahcall et al. 1982 
(No diffusion) 




19 




15.50 




13.6 


0.226 


4.3 




5.6 



TABLE IV. Central temperatures of 12 recently-computed solar models (no diffusion). 

Model Tc 
(106 K) 

Castellani, Degl'Innocenti, and Fiorentini (1993) 15.72 

Berthomieu et al. (1993) 15.55 

Turck-Chieze and Lopes (1993) 15.43 

Ahrens, Stix, and Thorn (1992) 15.65 

Bahcall and Pinsonneault (1992) 15.57 

Christensen-Dalsgaard (1992) 15.68 

Guenther et al. (1992) 15.53 

Guzik and Cox (1991) 15.40 

Proffitt and Michaud (1991) 15.71 

Sackmann et al. (1990) 15.43 

BahcaU and Ulrich (1988) 15.64 

Lebreton and Dappen (1988) 15.54 
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TABLE V. Characteristics of the energy profile. 



Model 



A (<? - <?peak> FWHM 

(kcV) (kcV) (keV) 



(kcV) (kcV) (kcV) (kcV) (keV) 



Bahcall-Pinsonncault 1992 
(Helium diffusion) 



Bahcall-Pinsonneault 1992 
(Helium diffusion) 



1.29 



Bahcall-Pinsonncault 1992 


1.28 


(No diffusion) 




Bahcall-Ulrich 1988 


1.28 


(No diffusion) 




Bahcall et al. 1982 


1.23 


(No diffusion) 





1.24 



Bahcall-Pinsonncault 1992 


1.23 


(No diffusion) 




Bahcall-Ulrich 1988 


1.23 


(No diffusion) 




Bahcall ct al. 1982 


1.18 


((No diffusion) 





Ground-state decay 
0.856 1.63 0.56 1.07 



0.858 



0.857 



0.821 



1.040 



1.62 0.55 1.07 



1.62 0.55 1.07 



1.56 0.53 1.03 



Excited-state decay 
1.048 0.97 0.24 0.73 



0.96 0.24 0.72 



1.039 0.95 0.24 0.71 

1.002 0.91 0.22 0.69 



2.59 
2.58 
2.58 
2.39 



10.0 
9.9 
9.9 

8.8 



0.23 
0.22 
0.22 
0.21 



2.86 


11.3 


0.04 


2.82 


11.1 


0.04 


2.82 


11.1 


0.04 


2.62 


9.9 


0.03 
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TABLE VI. The effect of vacuum neutrino oscillations on the energy shift, A(z/e) [see Eq. 
for electron neutrinos and on the energy shift , A(f^), for muon neutrinos is given as a function 
of vacuum mixing angles, 0v, and the phase-angle, (p, see Eq. (|111| ), evaluated at the peak of the 
energy spectrum at which neutrinos are created, gpeak = 862.27 keV. The last column gives the 
energy shift, A(e.sc.), that would be measured in an electron-scattering experiment. The survival 
probability, averaged over the neutrino profile, for an electron neutrino to remain an electron 
neutrino is given in column 3. 



sin^ 201 



2?!'(gpeak)/7r 



Survival {ve 



(keV) 



(keV) 



A(e.sc.) 
(keV) 



0.0 


0.0 


1.00 


3.5 


1.00 


4.0 


1.00 


4.5 


1.0 


4.6 


1.0 


4.7 


1.0 


4.8 


1.0 


4.9 


1.0 


5.0 


1.0 


5.1 


1.0 


5.2 


1.0 


5.3 


1.0 


5.4 


1.0 


5.5 


1.00 


6.0 


1.00 


7.0 


1.0 


7.2 


1.0 


7.4 


1.0 


7.6 


1.0 


7.8 


1.0 


8.0 


0.75 


3.5 


0.75 


4.0 


0.75 


4.5 


0.75 


5.0 


0.75 


6.0 


0.75 


7.0 


0.75 


8.0 



1.000 


1.286 




1.29 


0.495 


1.26 


1.31 


1.27 


0.9999 


1.29 




1.29 


0.507 


1.32 


1.26 


1.31 


0.352 


1.33 


1.26 


1.31 


0.212 


1.35 


1.27 


1.31 


0.100 


1.39 


1.27 


1.31 


0.027 




1.28 


1.30 


0.0002 




1.29 


1.29 


0.022 




1.29 


1.27 


0.091 




1.30 


1.26 


0.200 


1.22 


1.30 


1.26 


0.338 


1.24 


1.31 


1.26 


0.492 


1.25 


1.32 


1.26 


0.9997 


1.29 




1.29 


0.0004 




1.29 


1.29 


0.089 


... 


1.300 


1.25 


0.335 


1.22 


1.32 


1.25 


0.643 


1.25 


1.36 


1.26 


0.897 


1.27 


1.45 


1.27 


0.9995 


1.28 




1.28 


0.621 


1.27 


1.31 


1.28 


0.9999 


1.29 




1.29 


0.630 


1.30 


1.26 


1.30 


0.250 


1.29 


1.29 


1.29 


0.9998 


1.29 




1.29 


0.250 


1.29 


1.28 


1.29 


0.9996 


1.29 




1.29 



58 



